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Abstract 

We predict the existence of a new ferromagnetic shape memory alloy Ga2MnNi using density 
functional theory. The martensitic start temperature (Tm) is found to be approximately pro- 
portional to the stabilization energy of the martensitic phase (SEtot) f° r different shape memory 
alloys. Experimental studies performed to verify the theoretical results show that Ga2MnNi is 
ferromagnetic at room temperature and the Tm and Tq are 780 K and 330 K, respectively. Both 
from theory and experiment, the martensitic transition is found to be volume conserving that is 
indicative of shape memory behavior. 

PACS numbers: 71.15.Nc, 81.30.Kf, 75.50.Cc 
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I. INTRODUCTION 



Ni 2 MnGa exhibits both ferromagnetism and shape memory effect, and is a promising 
candidate for technological applications because of its high actuation frequency compared 
to conventional shape memory alloys.^ The unusually large strain caused by a moderate 
magnetic field (10% at 1 Tesla) 1 and the observation of giant magnetocaloric effect^ and 
large negative magnetoresistance^ in Ni 2 MnGa have started intense research activity in fer- 
romagnetic shape memory alloys (FSMA). However, a major drawback of Ni 2 MnGa is its 
brittleness. So, the present challenge in FSMA research lies in the search for new materials 
that have magneto-mechanical properties superior to Ni 2 MnGa, and preferably having high 
martensitic start temperature (Tm) and Curie temperature (Tc). In recent past, differ- 
ent groups have attempted to find FSMA materials with properties superior to Ni 2 MnGa. 
Takeuchi et al. have studied a range of compositions in the Ni-Mn-Ga phase diagram and 
found that the martensitic transition temperature decreases as the magnetization increases.- 
Although martensitic transition and inverse magnetocaloric effect have been reported re- 
cently in non-stoichiometric compositions of Ni-Mn-Sn, Ni-Mn-In, Ni-Co-Mn-In,-^ these 
systems have not emerged as viable alternatives to Ni 2 MnGa. 

Since related stoichiometric alloys like Ni 2 MnAl, Ni 2 MnIn, or Ni 2 MnSn do not exhibit 
martensitic transition,- it is apparent that Ga plays an important role in making Ni 2 MnGa a 
shape memory alloy. Substitution of Ga by In in Ni 2 MnGa decreases T^.- Thus, excess Ga 
may have a stabilizing effect on the martensitic phase. Zayak et al. theoretically studied the 
role of Ga 4p states in the stability of the martensitic phase of Ni 2 MnGa. 10 While considerable 
experimental work has been done on Ni and Mn excess Ni-Mn-Ga,— li^iMiii this is not the 
case for Ga excess Ni-Mn-Ga. Theoretical investigations so far have concentrated on the 
electronic structure in Ni 2 MnGa and related stoichiometric Heusler alloys and total energy 
calculations have been done to ascertain the stability of the martensitic phase.— ii&ilii&ii2i22 

Here, using spin-polarized, full potential ah initio density functional theory, we establish 
a method to estimate the martensitic structural transition temperature and predict possible 
existence of a new FSMA Ga 2 MnNi. A tetragonal martensitic phase with c/a= 0.83 is found 
to be lower in total energy (E to t) compared to the cubic austenitic phase. The martensitic 
phase total energy is lower by 55 meV/atom (=5E to t, i.e. the stabilization energy, which is 
the difference of E tot between the austenitic and martensitic phases). This value is larger than 



2 



other FSMA materials studied by us.— Based on our data and those in Refs. |23| and \24 , 
we show that SE to t is approximately proportional to the martensitic transition temperature 
(Tm). E tot for the ferromagnetic state is lower than the paramagnetic state, showing that 
Ga2MnNi is ferromagnetic. Inspired by the theoretical prediction, Ga2MnNi has been 
prepared, and it indeed exhibits a thermoelastic martensitic transition with Tm= 780 K, 
which is highest reported so far in the Ni-Mn-Ga family. The Curie temperature (Tc) 
is 330 K. X-ray diffraction (XRD) shows that the structure is monoclinic with b~7xa, 
indicating the existence of modulation and hence, the possibility of magnetic field induced 
strain. Both from experiment and theory, we find that there is hardly any unit cell volume 
change across the martensitic transition, and this is strongly indicative of shape memory 
behavior. 25 



II. METHODOLOGY 

The ab-initio relativistic spin-polarized full-potential-linearized-augmented-plane-wave 
(FPLAPW) method calculations were performed using WIEN97 code^ with the general- 
ized gradient approximation for exchange correlation. An energy cut-off for the plane wave 
expansion of 16 Ry is used {RMTK max = 9). The cut-off for charge density is G max = 14. The 
maximum I (l max ) for the radial expansion is 10, and for the non-spherical part: l m ax,ns = ^- 
The muffin-tin radii are Ni: 2.2488, Mn: 2.3999, and Ga: 2.2488 a.u. The number of k 
points for self-consistent field cycles in the irreducible Brillouin zone is 413 and 1063 in the 
austenitic and martensitic phase, respectively. E tot consists of the total kinetic, potential 
and exchange correlation energies of a periodic solid. 27 The convergence criterion for the 
total energy E tot is 0.1 mRy, which implies that accuracy of E tot is ±0.34 meV/atom. The 
charge convergence is set to 0.001. The tetrahedron method for the fc-space integration has 
been used. 

Polycrystalline ingots of Ga2MnNi were prepared by melting appropriate quantities of 
the constituent metals of 99.99% purity in an arc furnace under argon atmosphere and 
subsequently annealed in sealed quartz ampule wrapped in Mo foil at 873 K for 12 days, 
then at 723 K for 1 day and finally slowly cooled to room temperature. The differential 
scanning calorimetry (DSC) measurements were done by using TA Instruments MDSC model 
2910 at a scan rate of 10°/minute. Magnetization was performed using vibrating sample 
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magnetometers (VSM) from Oxford Instruments and Lake Shore Cryotronics, Inc. Powder 
XRD patterns were obtained using Cu Ka radiation with Rigaku XRD unit at a scan rate of 
2° /minute. Energy dispersive analysis of x-rays (ED AX) was done using scanning electron 
microscope with Oxford detector model with 2% accuracy by estimating the intensities of 
Ni, Mn and Ga Ka characteristic lines (5.9 to 9.2 keV) that are well separated and have 
small background by averaging over several measurements. 

III. RESULTS AND DISCUSSION 

A. Theoretical studies using FPLAPW method 

The positions of the atoms in the cubic austenitic phase of Ga 2 MnNi are determined from 
the E tot calculations in the L 2l cubic structure that consists of four inter penetrating f.c.c. 
lattices at (0.25,0.25,0.25), (0.75,0.75,0.75), (0.5,0.5,0.5) and (0,0,0) (Fig. la, b). The 
first two positions are equivalent (8/), whereas the other two are 4a and 46, respectively. 
In our notation, GaGaMnNi means that the two Ga atoms occupy (0.25,0.25,0.25) and 
(0.75,0.75,0.75) i.e. the 8/ positions, while Mn and Ni are at (0.5,0.5,0.5) and (0,0,0), re- 
spectively. Similarly, GaNiGaMn means that Ga atoms occupy inequivalent (0.25, 0.25, 0.25) 
and (0.5,0.5,0.5) positions, while Ni and Mn atoms are at (0.75,0.75,0.75) and (0,0,0), 
respectively. E tot has been calculated as a function of lattice constant (a) for all the differ- 
ent possible Ga positions (GaGaMnNi, GaGaNiMn, NiMnGaGa, NiGaMnGa, GaNiGaMn, 
GaMnGaNi), where the two Ga atoms occupy either symmetry equivalent or inequivalent 
points. Etot values for the inequivalent Ga structures (NiGaMnGa, GaNiGaMn and GaM- 
nGaNi, Fig. lb) are similar. The equivalent Ga structures (GaGaMnNi, GaGaNiMn and 
NiMnGaGa, Fig. la) are also very close to each other in energy. The data have been fit- 
ted using a least square minimization routine using the Murnaghan equation of state (solid 
lines, Fig. lc). The minimum E tot for the equivalent Ga structures (arrow) is lower by 113 
meV/atom compared to the inequivalent Ga structure (tick), unambiguously denoting the 
former to be the stable structure of Ga 2 MnNi in the austenitic phase. The E tot minimum 
(arrow) is at a= 11.285 a.u. (5.96 A) with the unit cell volume of 1437 a.u. 3 (Fig. lc). 
Furthermore, the formation energy of Ga 2 MnNi is calculated by £ , to t(Ga 2 MnNi)-2x£' tot (Ga)- 
£' to t(Mn)-i?t i(Ni). The formation energy turns out to be negative, comparable to Ni 2 MnGa, 
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indicating that the compound is stable. It should be noted that all the calculations shown 
in Fig. lc have been performed in the ferromagnetic state, since this is the stable magnetic 
phase (discussed later). 

The martensitic transition involves a structural transition from cubic to a lower symmetry 
phase with decreasing temperature. In order to study this phase transition in Ga2MnNi, our 
strategy is to calculate E to t as a function of a volume conserving tetragonal distortion by 
varying c/a. As c/a is increased from the cubic value of unity, E to t increases (Fig. lc). On 
the other hand, for c/a<l, E tot initially decreases and a minimum is obtained at c/a= 0.83 
(dashed tick). In the next step to reach the global E to t minimum in the martensitic phase, 
the unit cell volume is varied keeping c/a fixed, and the minimum is obtained at the unit 
cell volume of 1435.8 a.u. 3 with a— 12.004, c— 9.964 a.u. (Fig. 2a, dashed arrow). Thus, 
although there is a large change in lattice constants (+6.4% in a and -11.7% in c), there 
is almost no volume change between the austenitic and the martensitic phases. E tot has 
been calculated for Ga2MnNi in the paramagnetic state in the martensitic phase using the 
optimized lattice constants. It turns out to be 156 meV/atom higher than the ferromagnetic 
state. Thus, Ga2MnNi has a ferromagnetic ground state. The total spin magnetic moment 
of Ga2MnNi in austenitic (martensitic) phase is 3.04 (2.97) The local moments of 

Mn, Ni, and Ga in the austenitic (martensitic) phase are 3.03 (2.87), 0.06 (0.16) and -0.05 
(-0.05) hb, respectively. On the basis of the condition that a volume conserving martensitic 
transition is the necessary and sufficient condition for shape memory behavior— and that the 
ground state is ferromagnetic, we predict Ga2MnNi will behave as a ferromagnetic shape 
memory alloy. 

The martensitic phase being the lower temperature phase, E tot for the martensitic phase 
is lower than the austenitic phase by 55 meV/atom. Larger stabilization energy i.e. 5E to t 
would imply greater stability of the martensitic phase and enhanced Tm- From our ear- 
lier calculations, 5E tot (experimental Tm) is found to be 3.6 (210 K), 6.8 (270 K) and 39 
(434 K) meV/atom for Ni 2 MnGa, Mn 2 NiGa and Ni2.25Mno.75Ga, respectively^^^^ Here, 
we report similar calculations for N^Mnln and Ni2MnAl. The optimized lattice constants 
of Ni 2 MnAl and Ni 2 MnIn (5.79 and 6.06A, respectively) are in good agreement with ex- 
periment: 5.83 and 6.08A. 29 Although their off-stoichiometric compositions exhibit marten- 
sitic transition, it is well known that neither of these Heusler alloys undergo martensitic 
transition.- Interestingly for N^Mnln, 5E to t turns out to be almost zero (0.34 meV/atom) 
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within the theoretical accuracy limit, while 5E to t for Ni 2 MnAl is negative (-0.94 meV/atom). 
These values of SE to t indicate that the martensitic phase in Ni 2 MnAl and Ni 2 MnIn is not 
stable and so martensitic transition will not occur. This is in agreement with experimental 
data and earlier theoretical work.— ^ 

From the above data, a correlation emerges between SE tot and Tm- Conceptually, this 
is understandable since larger 5E to t implies higher stability of the martensitic phase at zero 
temperature. A first order transition to the austenitic phase would occur when with increas- 
ing temperature, the martensitic phase energy (defined by the energy minimum in Fig. 2a) 
would increase to reach the energy minimum for the austenitic phase. This means with 
increasing temperature, to undergo the martensitic transition, the energy of the martensitic 
phase has to overcome 5E tot and this would be directly related to k,BT M - A similar concept 
has been used in Ref. |2l|, where taking 5E to tO<:kBTM, the increase in Tm between Ni 2 MnGa 
and Ni 2 .25MnGa could be explained. This expression should be generally valid, and this 
indeed seems so for TiNi (45, 333), TiPd (95, 783) and TiPt (155, 1343). The numbers 
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in bracket indicate 5E to t and Tm in meV/atom and K, respectively, as taken from Ref. 
Similar trend is obtained for Ni excess Ni-Mn-Ga.— 

In Fig. 2b, T M versus SE tot for all the shape memory alloys discussed above are plotted; 
Tm is taken to be zero for Ni 2 MnIn and Ni 2 MnAl. It is highly significant that although 
the theoretical data are from three different groups^^i^ 1 ^ 1 ^ Q n two different types of 
shape memory alloys and the methods of calculation are different, an approximately linear 
relation between Tm and SE tot is evident. Thus the validity of the expression SEtotocksTM 
is established. A rather good straight line fit through the data for TiX (= Ni, Pd, Pt)^ is 
obtained (Fig. 2b). Since the Ni-Mn-X (X= Ga, In, Al) FSMA's are different from TiX, 
a separate straight line is fitted. The quality of the fit is similar to TiX; except for data 
around 200 K. This is possibly because of the existence of modulated structures this Tm 
range, which is not considered in theory. From the fitted line, Tm for Ga2MnNi is estimated 
to be about 570 K (filled circle), corresponding to its SE to t=55 meV (Fig. 2b). 

It is generally believed that Tm would increase with the valence electron per atom ratio 
(e/a). However, this relation is of limited applicability and breaks down in many cases: 
for example, Ni 2 MnGa, Ni 2 MnIn and Ni 2 MnAl all have the same e/a (=7.75), but only 
Ni 2 MnGa exhibits a martensitic transition. TiX (X=Ni, Pt, Pd) has the same e/a (=6.5), 
but their Tm is very different. In Ni-Mn-Ga-In, although e/a is same, Tm changes.- For 
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Ni2_a;Mni +a ;Ga between x=0.25 to 1, we find that as e/a decreases from 7.31 to 6.75, Tm 
increases from 37 to 270 K.— ^ For the alloys shown in Fig. 2b, the absence of any correlation 
between Tm and e/a is shown as an inset. In contrast, the present approach explains all 
the above observations. For example, SE tot decreases from 3.6 meV/atom to zero between 
Ni 2 MnGa and Ni 2 MnIn, which explains the decrease in Tm with In doping and the absence 
of a martensitic transition in N^Mnln. Higher SEtot i n M^NiGa rationalizes why its Tm is 
higher than N^MnGa, although its e/a (=6.75) is lower. Thus, the proportionality of Tm 
with SEtot is of more general validity, since it has a theoretical foundation that involves all 
electron ab-initio calculations, unlike the phenomenological relation between Tm and e/a. 
In fact, this approach to determine the transition temperature should be applicable to any 
first order structural transition. 

B. Experimental studies 

Differential scanning calorimetry on polycrystalline ingots of Ga2MnNi shows a clear 
signature of a first order martensitic transition with Tm= 780 K and austenitic start tem- 
perature (A s ) of 790 K (Fig. 3a). The experimental Tm is considerably higher than the 
theoretically predicted value, and a possible reason is discussed below. The latent heat of 
the transition turns out to be about 2.35 KJoule/mole, which is similar to that reported 
for Ni excess Ni-Mn-Ga, for example, Ni2.24Mno.75Ga.— The difference in the width of the 
heating and cooling thermograms could be related to the kinetics of the structural tran- 
sition. EDAX measurements from different regions of 30/ix30/i well as the back 
scattered image show that the specimen is homogeneous. The average composition turns 
out to be Ga1.9Mn1.08Ni1.02- In agreement with theory, the isothermal M — H curve at 
2.5 K shows that Ga2MnNi is indeed ferromagnetic (Fig. 3b). The hysteresis loop is not 
clearly observed because the coercive field is small (^25 mT). Such small coercive fields have 
been reported for other Ni-Mn-Ga alloys.— 1 ^ The saturation field is 1 T and the saturation 
moment is 1 //e/f.u. M(T) at low field gives 7^=330 K (arrow, Fig. 3c). This implies that 
the martensitic transition occurs in the paramagnetic state and expectedly M(T) shows no 
change across Tm- If should be noted that the saturation moment of 1 /xg/f.u. is less than 
the theoretically calculated moment of about 3/x^/f.u. The reasons for this disagreement 
could be that the actual sample has Mn excess, which might cause Mn clustering leading 
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to antiferromagnetic coupling between Mn atom pairs, as has been observed for other Mn 
excess systems.— Moreover, note that the theory does not consider the actual mono- 
clinic structure (discussed below) which might favor a different magnetic ground state with 
anti-parallel coupling between Mn atoms. 

The x-ray diffraction (XRD) pattern corresponding to the austenitic phase has been 
simulated by the Le Bail fitting procedure, and the structure is clearly cubic L 2l . The 
relative intensity of the (200) peak compared to the (111) peak (shown in an expanded 
scale in Fig. 3d) confirms that the Ga atoms occupy the equivalent 8/ position, in agree- 
ment with theory (Fig. 1). The experimental lattice constant (a aus =5.84 A) is close to 
the calculated value (5.96 A). However, the martensitic phase XRD pattern is more compli- 
cated than tetragonal and can be indexed by a monoclinic phase (P2/m space group) with 
a= 4.31, b= 29.51 and c= 5.55 A, and f3— 90.49. Since 6~ 7xa, a seven layer modulation may 
be expected, and such structures with monoclinic or orthorhombic symmetry that exhibit 
modulation has been reported for Ni-Mn-Ga.— Magnetic field induced strain has been ob- 
served in Ni-Mn-Ga for structures that exhibit modulation.- The c/a for this monoclinic cell 
(that can be compared to the theoretical c/a= 0.83 for the tetragonal structure) is obtained 
by c/a= 5.55/(4.31x \/2)= 0.91. Thus, the agreement between experimental and theoretical 
c/a is reasonable, considering that a simplified structure is used in theory. 

However, the most important point is that the experimental unit cell volume of the 
martensitic phase is within 1% of that of a comparable austenitic cell given by 7xa^ us /2. 
This shows that the unit cell volume hardly changes between the two phases, which is a nec- 
essary condition for a shape memory alloy. Thus, a unit cell volume conserving martensitic 
transition with small width of hysteresis (Fig. 3a) and presence of modulation indicate that 
Ga2MnNi is indeed a FSMA material. 

IV. CONCLUSION 

The modulated martensitic structure of Ni-Mn-Ga is complicated and a controversy exists 
even about the structure of the well studied Ni 2 MnGa. 34 Atomic positions have not yet been 
determined for the monoclinic structure. Under such circumstances, our work is important 
because it shows that a new FSMA material can be predicted by computing the energy 
cost of formation of the martensitic phase in a simpler tetragonal structure. The present 
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work demonstrates that a new FSMA material can be predicted by determining the energy 
stability of a tetragonal martensitic phase with respect to the cubic austenitic phase. This 
approach is successful because, although the modulated phase involves a large unit cell, the 
atoms are generally displaced only by a small amount from their positions compared to the 
tetragonal structure.— Since the tetragonal structure is not computationally demanding, 
precise calculations can be performed for lattice constant optimization in the lowest energy 
magnetic state.— i2L22 Thus, the total energy difference can be determined with sufficient 
accuracy and thus Tm estimated. However, difference in Tm between experiment and theory 
could occur, as in this case, possibly because the latter does not consider the actual structure. 
In this context, it is to be noted (Fig. 2) that a subtle change in 5E tot can substantially alter 
the Tm value. Theory thus provides an important starting point for the experimentalists, 
and experimental inputs can be used to further refine the theory. A direct proof of the 
FSMA behavior is the movement of twins with magnetic field and the actuation behavior. 
So, further work on the magneto-mechanical behavior of Ga2MnNi is in progress. Prediction 
of new materials in the quest for better properties is the need of the hour in FSMA research 
and the present work aims towards that. 
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